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We report the discovery of a metamagnetic phase transition in a polar antiferromagnet NisTeOe that occurs 
at 52 T. The new phase transition accompanies a colossal magnetoelectric effect, with a magnetic-field-induced 
polarization change of 0.3 pC/cm^, a value that is 4 times larger than for the spin-flop transition at 9 T in the 
same material, and also comparable to the largest magnetically-induced polarization changes observed to date. 
Via density-functional calculations we construct a full microscopic model that describes the data. We model the 
spin structures in all fields and clarify the physics behind the 52 T transition. The high-held transition involves 
a competition between multiple different exchange interactions which drives the polarization change through 
the exchange-striction mechanism. The resultant spin structure is rather counter-intuitive and complex, thus 
providing new insights on design principles for materials with strong magnetoelectric coupling. 


Magnetoelectric (ME) multiferroics have been extensively 
studied recently to understand the mechanisms responsible 
for cross-coupling between magnetism and ferroelectricity, 
which is at the heart of their promise for application in mul¬ 
tifunctional devices w. In this class of materials, at least 
three mechanisms are known to induce ferroelectric polariza¬ 
tion (P) upon magnetic order: (1) the spin current or inverse 
Dzyaloshinskii-Moriya (DM) interaction in a spin-cycloidal 
structure which is mediated by anti-symmetric exchange ||7E 
|9], (2) the symmetric exchange-striction mechanism between 
parallel or anti-parallel alignment of spins ifTOl . and (3) the 
hybridization between metal d- and ligand p-orbitals that is 
modulated by spin direction ca. A majority of ME cou¬ 
plings that have been studied to date involve mechanisms (1) 
and (3). However, the symmetric exchange mechanism (2) 
can, in principle, also lead to large ME effects. 

Another route to a large ME effect is to consider mag¬ 
netic systems that have a polar structure. This condition 
meets the prerequisites for the ME effect, i.e., coexistence of 
broken spatial-inversion-symmetry and time-reversal symme¬ 
try. Often, these systems exhibit a non-polar to polar struc¬ 
tural transition at high temperatures and magnetic ordering at 
lower temperatures. A well-known example is BiEeOa, with 
ferroelectric and antiferromagnetic transition temperatures at 
1100 K and 653 K respectively ifT^ . However, its ME effect 
is small compared to those of spin-driven ME materials. 

Recently, a ME effect has been observed in the corundum- 
related compound NiaTeOe (NTO). It crystallizes in a po¬ 
lar R3 space-group with three Ni^+ ions (3d®, S' = 1) and 
a non-magnetic Te ion arranged along the c-axis in a unit 
cell (Eig. [TJa)) to form a linear chain with broken spatial- 
inversion symmetry. This material is not ferroelectric but py¬ 
roelectric, i.e., the electric polarization cannot be switched by 
an external electric field. However, in addition to its non- 
switchable electric polarization component due to the polar 
structure, it also shows a large magnetically-induced polariza¬ 
tion. In zero magnetic field (H), below the Neel temperature 
Tn = 52 K, NTO has a collinear antiferromagnetic (AEM) or¬ 



FIG. 1. (a) Crystal structure of NisTeOg. Only Ni ions are shown. 
(b)-(e) Schematic spin structure along a c-axis chain at different mag¬ 
netic fields applied along the c-axis. The spins in the buckled hon¬ 
eycomb planes are aligned ferromagnetically in zero magnetic field 
ini. Ji stands for the five nearest-neighbor exchange constants con¬ 
sidered in the text. Dotted boxes indicate the magnetic unit cell along 
the c-axis, which is doubled for (b)-(d) but not for (e). 


der, tttiii in the rhombohedral unit cell with spins aligned 
along the c-axis ini. The spatial inversion symmetry is bro¬ 
ken and the electric polarization is pointing along the c-axis, 
as in Ca3(Co,Mn)06 ifTOlflTll . It was found that NTO under¬ 
goes a second-order spin-flop (SE) transition at a critical field 
Hci 9 T, which accompanies a large ME effect ifTSll . Here, 
symmetric exchange-striction at the SE transition distorts the 
polar crystal structure to modify the electric polarization. The 
ME coefficient (a = ^) is as high as 1300 ps/m at Hd, with¬ 
out any magnetic hysteresis. However, a phenomenological 
model with only two magnetic sublattices was implemented 
to describe the polarization change at the SE transition flSl . 
whereas a full description of NTO requires a model with three 
different Ni spins and all possible exchange interactions be¬ 
tween them lfT3]l . 

In this Letter, we present a new phase transition in NTO 
discovered by high magnetic held study up to 92 T. The high 
held transition is accompanied by a colossal ME effect that is 
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FIG. 2. (a), (b) Magnetization and (c), (d) change of electric polar¬ 
ization (AP) of NisTeOe for magnetic fields applied along different 
crystalline axes as indicated, (e) c-axis magnetostriction as a func¬ 
tion of magnetic field applied along the c-axis. Solid lines are ex¬ 
perimental data taken under isothermal condition and open squares 
are from model calculations based on parameters shown in Table 
(f) AP measured up to 92 T for H || ab and H |j c. (g), (h) The 
magnetic field-dependence of M and AP along the c-axis in P || c 
configuration obtained from the model calculations described in text. 
Dashed lines in (a), (b) are guides for the eyes, (a), (b) magnetiza¬ 
tion and (e) magnetostriction curves are shifted for clarity, (e) Jump 
in magnetostriction at Hc 2 shows a slight difference in magnitude 
between different types of magnet used, possibly due to the variation 
in magnetic field sweep rate, (f) Lines and open symbols denote data 
taken by using a capacitor-bank-driven and a hybrid pulse-magnet, 
respectively. A single domain sample was used for (c), (d) while a 
multi-domain (as-grown) sample was used for (f). 


comparable to largest field-induced polarization changes re¬ 
ported to date GSiini. We corroborate our experimental re¬ 
sults with a microscopic model, using model parameters ex¬ 
tracted from extensive density-functional theory (DFT) calcu¬ 
lations, thereby predicting the magnetic structure at all mag¬ 
netic fields. Multiple exchange interactions are found to con¬ 
tribute to the field-dependence of electric polarization via an 
exchange-striction mechanism. 

Fig. I^a) shows the magnetic field-dependence of the mag¬ 
netization M along the c-axis of NTO up to 92 T. At 4 K 
a sharp SF transition is evident at Hci 9 T, then M in¬ 
creases quasi-linearly up to 20 T. The linear extrapolation of 
the M(H) data between 9 and 20 T has a zero intercept at 
H = 0, consistent with a SF transition. When H is further in¬ 


creased above 30 T, the slope of the M(H) curve decreases 
slightly. At Hc 2 = 52 T, a small jump in M is observed. 
With further increasing H, M increases linearly up to 92 T. 
The value of M at 92 T is 4.8 /is per formula unit (f.u.), 
which is still smaller than the expected saturation magneti¬ 
zation Ms = 6 fiB/f-Vi. for three Ni^+ ions with S' = 1 (assum¬ 
ing gyromagnetic ratio of 2). By assuming that there are no 
other magnetic field-induced transitions, linear extrapolation 
of the M(H) curve to the expected Ms value gives a satu¬ 
ration magnetic field Hs of 120±10 T. The in-plane M(H) 
data (Fig. |^b)) below 30 K show a monotonic increase un¬ 
der magnetic field. At 30 K, the M(H) curve shows a cusp 
at iT = 60 T which decreases upon further warming. The 
most striking feature of the high-held transition in NTO is 
the colossal change of P at Hc 2 and reversal of the AP{H) 
slope at higher fields. Fig. |^c) shows the change of c-axis 
electric polarization AP = PiH) - P(H=Q) as a function 
of magnetic field applied along the c-axis. In this configu¬ 
ration, AP(H) curve shows a step at Plci as previously re¬ 
ported M- When the magnetic field is further increased at 
4 K, P slightly increases, and then shows a sudden drop at 
50 T, close to Hc 2 , with a minimum at 53 T. The overall AP 
at iTc 2 reaches 0.3 pC/cm^ at 4 K. Counter-intuitively, AP at 
Plc 2 is ~ 10 times larger than that at Hd, whereas the change 
of M at Hc 2 is almost two times smaller than that at Hd- 
When the magnetic field is further increased above 53 T, P 
increases linearly up to 65 T. In contrast, the c-axis polariza¬ 
tion measured under in-plane magnetic field does not show 
any sharp jump, but only a smooth reversal in slope which 
shifts to lower magnetic field as the temperature is increased 
(Fig.l^d)), concurrent with a cusp in the in-plane M{H) curve 
(Fig. rab)). The lattice also responds sensitively to the mag¬ 
netic field at these transitions, as shown in the c-axis magne- 
tostriciton AL/L measurements (Fig.j^e)). 

We further explored AP(iT) up to 92 T, with a different, 
multi-domain sample (Fig. |^f)). When the magnetic field is 
applied along the c-axis, P increase linearly above Hc 2 , up 
to 92 T. A linear extrapolation of AP{H) curve above Hc 2 
gives AP(iT) = 0 at 120±5 T, consistent with the expected 
saturation magnetic field from the M{H) curve. This im¬ 
plies that the magnetically-induced electric polarization is no 
longer active when the system is in the fully saturated phase. 
At elevated temperatures, the sharp changes of P at 11^2 are 
still observed, although its magnitude and transition field de¬ 
creases up to 30 K, above which the sharp drop is replaced 
with a slope change in the AP(H) curve. When the magnetic 
field is applied along the a6-plane, the AP(H) curve shows 
only a sharp reversal of slope at 70 T at 4 K. Above 70 T, 
two AP(H) curves measured in different configurations co¬ 
incide with each other, suggesting an isotropic magnetic be¬ 
havior above this field. 

The field-induced AP value of 0.3 pC/cm^ (Fig. |^c)) and 
the high ME coefficient value ESI at Hc 2 in NTO are among 
the largest observed in magnetoelectric systems EUElEll. 
All spin-driven ME materials exhibit a change of P when 
a field-induced phase transition or spin-reorientation occurs 
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FIG. 3. H-T phase diagram of NisTeOe with magnetic field 
applied along the c-axis determined by peak positions in dM/dH, 
a = dP/dH, and specific heat (see Supplement 1181 1 data. Dashed 
line is a guide for the eyes. 


ESIISI. However, for most of them, AP(H) is typically 
less than 0.01 /iC/cm^. Thus, NTO is a prototypical exam¬ 
ple where the polar symmetry and the additional polarization 
coupled to the magnetically ordered state give rise to a large 
ME coupling. 

The field-dependence of AP is similar to the case of 
BiFeOa along certain directions in that the AP(H) curve 
shows a step-like feature and a slope change EolEI]. How¬ 
ever, in the case of BiFeOa the field-dependent behavior is 
due to a phase transition from a spin-cycloid to a canted-AFM 
state, whereas we do not find any evidence of spin-cycloid or 
spiral structure in NTO. In addition, the spin-spiral state in 
NTO is only allowed in the a&-plane by symmetry, where it 
cannot contribute to c-axis polarization. 

We note that AP displays no electric field-dependence at 
either magnetic transitions ca, suggesting that NTO is not 
a type-II multiferroic where magnetic order induces ferro- 
electricity 0. Rather, the magnetic order modifies the exist¬ 
ing electric polarization associated with the polar space-group 
which is established already at very high temperature (a pos¬ 
sible ferroelectric Tc of ~ 1000 K was reported E^l l. 

Using our experimental results, we construct an H-T phase 
diagram for NTO in magnetic field along the c-axis, shown in 
Fig.[^ We observe three ordered phases below AFM, SF, 
and metamagnetic phase (MM) which are separated by phase 
boundaries, determined by M(H) and AP{H) curves, that 
are almost vertical at low temperatures. The high-temperature 
phase boundary between the paramagnetic (PM) and the or¬ 
dered phase (dotted line in Fig. extrapolates linearly to 
120 T at r = 0, again, consistent with the extrapolation of 
the M(H) and AP(H) curves to their expected saturation and 
zero values, respectively. 

Turning now to the theoretical modeling of this material, 
we note that the phenomenological description of magnetism 
in NTO lITSl is applicable in the vicinity of the low-field SF 


transition, but may not be accurate away from it. In order to 
study magnetic transitions in the whole magnetic field range, 
we use a simplified microscopic model with the Hamiltonian 


• 4 + E [-K2AS^,.f -H-Si 


, ( 1 ) 


where we model the Ni spins Si (S = 1) classically, and i,j 
run over all Ni sites in the lattice. The exchange constants be¬ 
tween Ni spins are taking particular values Ji... for 
the bonds (ij) indicated in Fig. We neglected anisotropic 
exchanges as well as next-nearest-neighbor Heisenberg ex¬ 
changes, since the second transition appears in the model 
without them. The terms with K 2 and H model single-ion 
easy-c-axis anisotropy and the coupling to a uniform magnetic 
field, respectively. 

The exchange constants depend on the ionic coordinates, 
and as a result, the ions shift in response to spin reorienta¬ 
tions in such a way as to strengthen the exchanges that favor 
the existing spin arrangement. These shifts of charged ions in 
a polar structure result in a change of P, which, assuming 
Heisenberg exchange-striction dominates, can be expressed 
(neglecting higher-order terms containing {Si ■ Sj)'^) as 

AP, = Y(^nSn-S'„, ( 2 ) 


where Sn and S!^ are the spins connected by the exchange 
interaction Jn, and the are exchange-striction parameters. 
We use Eq. Q to model the dependence of the polarization P 
on the spin configuration. 

The coefficients = an,ion + ctn.eU with the two terms 
describing the polarization contributions due to ion shifts 
and deformations of electronic orbitals, respectively, are cal¬ 
culated using DFT ifT^ . Similarly, the magnetostriction 
- the change of the sample size under an applied mag¬ 
netic field - is described by the equation, analogous to Q: 

The determination of exchange constants J„ is a difficult 
problem. The values calculated previously using DFT ||23]| 
give a non-collinear ground state when the energy is mini¬ 
mized within a magnetic unit cell at realistic values of K 2 . 
We have found that the exchange constants estimated using 
the PBEO hybrid functional approximation to DFT ll^ give 
the correct ground state and reproduce the second transition. 
We have then adjusted these constants to better fit the exper¬ 
imental M(H) data measured along the c-axis. The result¬ 
ing Jn parameters are summarized in Table along with the 
exchange-striction constants and magnetostriction param¬ 
eters Xn, calculated using DFT + U, as described in the Sup¬ 
plement im. 

With these model parameters in hand, we computed the 
changes in the spin arrangement under magnetic field applied 
along the c-axis; the results are shown in Fig. [TJb-e). The 
resulting magnetization and polarization curves, presented in 
Fig.|^g,h), are in qualitative agreement with experiments. We 
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TABLE I. Exchange (J„), easy-axis anisotropy constants (K 2 ), 
exchange-striction parameters for electronic (a„,ei) and ionic con¬ 
tributions (Qn.ion) to the electric polarization, and magnetostriction 
parameters (A„) estimated using DFT calculations and adjusted to 
ensure an AFM ground state. 


n 

1 

2 

3 

4 

5 

K 2 (meV) 


-0.6 

-3.1 

2.2 

4.2 

1.0 



-0.6 

-3.1 

2.2 

4.2 

0.69 


Q^n,el 

0.25 

3.3 

-0.1 

-3.0 

-0.8 


Q^n.ion 

-2.4 

-1.6 

-2.0 

11.2 

6.2 


An X 10'= 

2.1 

-1.6 

-4.3 

-8.9 

-15.3 


J™“(meV) 

-1.13 

-2.97 

0.79 

2.76 

0.32 

0.05 

J™“-"‘‘j(meV) 

-0.69 

-3.63 

0.76 

3.26 

0.65 

0.1 


checked that the calculated transition sequence and spin struc¬ 
tures did not change significantly with tuning of the exchange 
constants, suggesting an absence of competing phases. Our 
confidence in our correct identification of the phase transi¬ 
tions is further reinforced by the agreement of the calculated 
and measured magnetization and polarization curves. 

In the textbook SF transition for a two-sublattice antifer- 
romagnet, the exchange favors the collinear state ll25l . Sur¬ 
prisingly, in NTO the situation is the opposite - the AFM ex¬ 
change J 5 actually favors the canted state. Ferromagnetic ex¬ 
change J 2 favors the collinear state, in which the spins of Nil 
and Ni2 are parallel to each other, as shown in Fig.|^ The evo¬ 
lution of energy contributions from different exchanges un¬ 
der the applied magnetic field is illustrated in the Supplement 
ifTSl . In the canted state, above Hd, the spins of Nil and Ni2 
are no longer parallel thus losing energy on J 2 , but this cant¬ 
ing allows a gain in energy from other exchange interactions, 
while gaining Zeeman energy from the canting of the Ni2 and 
Ni3 spins along the magnetic field, as shown in Fig.|^c). 

As the magnetic field is increased further above Hd, the 
spin of the Ni2 rotates. At high-fields the c-component of the 
spins is pinned by the field, while the a 6 -plane components 
are chosen to minimize the energy (Eq. ([T]i). This is similar to 
the energy of the collinear state, except that the a 6 -component 
is not constrained to have the unit length. That is why the state 

that results above Hd, with six spins pointing, —->^ 0 ^, 

differs from the zero-field state, ttiiit- Above Hc 2 , the spin 
of the Nil cants further along the magnetic field, and simul¬ 
taneously the a 6 -plane components of the spins in half of the 
magnetic unit cell reverse in order to gain energy on the anti¬ 
ferromagnetic J 5 exchange, acting between Ni 1 andNiS spins 
from the neighboring crystallographic unit cells. At the same 
time the energy contributions from all the other exchanges in¬ 
crease, as evidenced by the total energy calculation ifTSl . A 
reversal of the a5-plane components of the spins in every sec¬ 
ond crystallographic unit cell leads to a large change of the 
magnetically-induced electric polarization and strains, as seen 
inFig.|gc,e ,h). We can see a restoration of translational sym¬ 
metry along the c-axis that was broken by AFM ordering. 

While the magnetic single-ion anisotropy plays an impor¬ 
tant role for the SF transition, the second transition at Hc 2 is 


controlled by the magnetic exchanges and the external mag¬ 
netic field, and is not sensitive to the single-ion anisotropy, 
thus it is not a classical SF transition. Our model predicts the 
transition at iJc 2 for both H W c and H || ab. Experimen¬ 
tally, however, no sharp transitions are observed in iT || ab, 
but there exists a cusp in the P(H) at around 70 T at T = 4 K, 
suggesting that the transition is of the second-order. As seen 
in Eig. |ga,b ,f), the magnetization curves and c-axis electric 
polarization with magnetic field applied different directions 
nearly coincide above the transition, implying that the mag¬ 
netic states for H \\ c and H \\ ab above the transition are 
similar. This discrepancy between the theory and experi¬ 
ment is partly due to the presently neglected DM interac¬ 
tions and symmetric anisotropic exchanges. In particular, our 
DET calculations indicate the presence of strong Ising-type 
anisotropies, which will modify the phase diagram for the c- 
axis ordered structures much more than the a&-plane ones. 
A detailed description of this aspect requires the introduction 
of additional parameters in our model, and additional experi¬ 
ments are being conducted in order to determine the parame¬ 
ters reliably. This work will be reported elsewhere. 

In summary, a high-held study of NTO reveals the pres¬ 
ence of a second spin reorientation transition well above the 
SE transition at 9 T. The high-held transition is hrst-order at 
52 T for H || c, and is second-order at 70 T for H |j ab at 
base temperature. Our theoretical analysis suggests that this 
high-held transition is governed by the competition between 
the Zeeman and exchange energies, and entails the reversal 
of the a&-plane component of half of the spins, leading via 
the Heisenberg exchange-striction to a change of electric po¬ 
larization that is among the largest observed to date. This 
spin reorientation results in restoration of translational sym¬ 
metry along the c-axis that was broken by AEM ordering. 
Our calculations also identify particular exchange interac¬ 
tions that are responsible for the stabilization of the magnetic 
phases in NTO. Eurthermore, compared to its isostructural 
compounds (MnsTeOe and CoaTeOe), which show multifer- 
roic behavior via antisymmetric exchange Il26l428ll . NTO ex¬ 
hibits an exchange-striction-driven polarization response that 
is almost two orders of magnitude larger, thereby demonstrat¬ 
ing a unique behavior in this class of materials. We propose 
that this scenario can form the basis of a new strategy for the 
design of materials with large ME effects. 
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